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Abstract 

> . 

CNJ , We use a solvable model to perform modified Dyson mapping and reveal the 

unphysical-state effects in the original Random Phase Approximation (RPA). 
We then propose a method to introduce the RPA and improve it based on the 
Boson description. 
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! I. Introduction 



The RPA is often used to describe collective excitations. In nuclear physics it is widely employed to 
define surface oscillation and pair oscillation fl. Recently the RPA has also been used to study mesons in 
^ ' terms of the quark-antiquark degrees of freedom^, [|. 

The RPA builds bosonic excitations with fermions so it includes unphysical state effects. In view of its 
^ ' important role in practice, it is helpful to explore convenient approaches for improvement. In this paper 
• <-h . we will develop such a method by studying a simple solvable model with the help of the modified Dyson 
' mapping^, We will introduce the RPA in a different way and explain it as the zero order approximation 
£3 . of a general Boson description. Then only very low orders of the correction terms will be included so as to 
preserve the simplicity of the method. It will be seen that our method are both quite powerful. Although 
the model is very simple our method can be, in principle, generalized approximately. 



II. Model and Boson Description for Particle- Antiparticle Pair States 

Consider a simple system of identical Fermions and their anti-particles and suppose that each particle 
has only two independent states. The Hamiltonian is assumed to be: 

H = 6 £ ( b l b » + d t d ^ + R(bUUbli4 + di6_id_i6i) , 

where 6t|0) and dfjO) stand for the particle states and anti-particle states respectively. The state vector of 
the ground state of the system is 

|*o> = Xl \0) + x 2 b\dl 1 bl 1 d\\0) , 

with 

where Eq is the ground state energy: 

So = 2e- v /(2e) 2 + i? 2 . 



The particle-antiparticle pair states bJdljO) and b^djlO) are also the exact eigen states of H belonging to 
the eigenvalue 2e. Namely, 

HbldL^O) = 2e&t<£ 1 |0), 
Hb^dllO) = 2e6Li4|0>- 
On the other hand the RPA excitation of a particle-antiparticle pair takes the following form 

QLi = Xd\bl 1 +Yb 1 d^ 1 , 
where X,Y are real numbers and to be found from the RPA equation: 

f [H,QI]kl>rpaQI 

I [Q,*,Qi] « <W 



The solution is : 

(x 2 - r 2 = 1) . 



J w MP = V(2e) 2 - fl 2 

wjj4P is the excited energy with respect to the RPA ground state which is determined by the condition: 

Q^ PA ) = 0. 



One gets from this 
with 



\^ PA ) = x^o) + x 2 b\dl 1 bl 1 d\\Q) • 



£-4--lw- 1 <* + *- 1 >- 



It can be seen clearly that the RPA worsens when ^ — > 1 and breaks down when ^ < 1 . 

Since the RPA operators , are assigned to satisfy the standard Boson commutation rule, it 
should be convenient to analyse the RPA within an exact Boson description. For this purpose, it is enough 
to consider the Boson image of the invariant subspace Vf of H, which is spanned by the states: 

\o) A<f-i\o) Ad-M4\o) ■ 

We now employ the modified Dyson mapping to set up a 1-1 correspondence between the fermion state 
subspace Vf and a Boson state subspace Vb- The latter is spanned by the four Boson states 

|0), Bt|0), BlilO), SjBlilO), 

where B^ are Boson creation operators, and are the corresponding annihilation operators. They of course 
satisfy the standard Boson commutation rule. These four Boson states arc to be described as the Boson 
images of the above four Fermion states. Therefore, |0) has the same quantum numbers as |0) and satisfies 

B„|0) = 0, (0|0) = 1, 

and B{\0) and Bl^O) carry the quantum numbers of &{<ili|0) and &L 1 d{|0) respectively. The correspondence 
can be expressed as 

$ = r B c/*, 

where 

U = (0|e B i t,, - l6l+B -i' ,l6 - 1 |0) ) 
t/t = (oig&I^Bx+fcl^Iu-iio). 



is the projection operator to the subspace V B . Namely: 

f $ when $ e Vb 
B ~ { when $ _L V B . 

According to this correspondence, one can easily find the Boson image of the Hamiltonian and write it as 

H B = r B {2e]T B^ + R(bIb^_ 1 + B^B^Tg ■ 

This Boson Hamiltonian and the Vb construct an exact Boson description for the particle-antiparticle pair 
states. Each Boson state of the subspace V B is a physical one, and a Boson state with components outside 
Vb is partially or totally unphysical. In this description, the 1-Boson eigenstates of Hb are B^O) which 
correspond to the 1-pair Fermion states, and the ground state is: 

|$o) =xi|0)+x 2 B 1 t B^ 1 |0). 

III. An Improvement on RPA Based on the Boson description 

The excitation elements of RPA are expressed in terms of Fermion operators and assumed to satisfy the 
standard Boson commutation rule. This cause some unphysical state effects which can clearly be seen from 
a Boson description. The exact Boson Hamiltonian H B can be written as 

H B = r B (wo + lurpa ^ A^A^Tb = T B H A r B , 

where 

lu = -2e + v / (2e) 2 - R 2 , 
4 = XBl + YB_, , 

H a =uj + ujrpa ^2 A^Aft . 

If we treat the operator Ha as a true Hamiltonian and solve its eigen equation in the whole Boson state 
space than we get the RPA results. Such a kind of eigenstates of H A of course contain many unphysical 
components. For a further study we now rely on the Boson description to introduce the RPA in a new point 
of view. Namely, we regard the following as the Boson Hamiltonian of RPA : 

H§ PA = T b (luo + lurpa B l B n) V B ■ 

Next, we will explain the RPA as the zeroth order approximation of our exact Boson description. For this 
purpose we make the transformation : 

f Bt=^t+aB_„ (/3 2 -a 2 = l), 
1 Bl = /?Bt-aB_„ 

and write H B as 

H B = T B {E' B t B » + R'( B t B -i + S_iBi)}r s , 

with 

E' = 4ea 2 - 2a/3R , 
J = 2e(a 2 + P 2 ) - 2a(3R, 
R' = (a 2 + f3 2 )R-4af3e. 



Introducing S and 7 

S = e -"i( B l B -i- B -i B i) ) 
e' 1 = a + (3, 
e - 7 = /3-Q!, 

one has 

St = 5Bt 5t , 

H B = T B S{E' W ]T B t B n + R'( B l B -i + 5_iBi)}5tr B . 

If a and are set to be X and Y respectively then one has 

Hb = FbSrpa(u + uj R pa B 1 B ^)^rpa^b ■ 

In other word, the RPA Hamiltonan H% PA can be found from the expression for Hb by replacing the 
operator S with 1 and then determining a and (3 with the help of the condition R' — 0. Our method to 
improve the RPA is to expand S in powers of 7 and keep only a few terms to get an approximate Hamiltonian: 

H B » T B {E a +ujY. B t B ^ + B ( B t B -i +B- 1 B 1 )}r B . 

Then the parameters are determined by the condition R = 0, and the following formulae: 

E~ = E' n + 2 7 i?' + 2 7 2 o)' + • • • , 

uj = J + 2jR' + 2~f 2 cu' H , 

R = R' + 2"/lu' + 2 7 2 i?' + • • • . 

We expect that our simple scheme (MRPA) will be useful in improving the behaviour of the PRA, especially 
in the region where the RPA begins to break down. Numerical results have confirmed our expectation. We 
present here results to second order of 7. We have also calculated the excitation energy to the fourth order 
and found it is a little better. This shows our method is stable and self-consistent. 

The figure is the one Boson excitation energy versus The energy unit is R. In order to show more 
clearly the improvement when — > 1, the abscissa is assumed to stand for ^ — 1 and is scaled logrithmically. 
When ^ — > 00, the results of the RPA and our method are close to the exact solution. When ^ = 1, both 
the RPA and our method give zero excitation energy. In the region where ^ is larger than 1 but not too 
large, the one Boson excitation energy of RPA deviates from the exact solution significantly and approaches 
zero rapidly while our method gives a much better description. The trend our method approaches zero is 
very slow which is an important improvement. 
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